Combination of a magnetic Feshbach resonance 
and an optical bound-to-bound transition 
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We use laser light near resonant with an optical bound-to-bound transition to shift the magnetic 
field at which a Feshbach resonance occurs. We operate in a regime of large detuning and large 
laser intensity. This reduces the light-induced atom-loss rate by one order of magnitude compared 
to our previous experiments [D. M. Bauer et al. Nature Phys. 5, 339 (2009)]. The experiments are 
performed in an optical lattice and include high-resolution spectroscopy of excited molecular states, 
reported here. In addition, we give a detailed account of a theoretical model that describes our 
experimental data. 

PACS numbers: 34.50.Cx, 33.20.Kf, 33.40.-|-f, 03.75.Hh 



I. INTRODUCTION 

Many properties of ultracold gases are determined by 
the interparticle interaction which is characterized by 
the s-wave scattering length a. This makes it desir- 
able to tune this parameter. A much-used method for 
this purpose is a magnetic Feshbach resonance [H, [3, H] . 
An alternative method is a photoassociation resonance, 
which is sometimes also called optical Feshbach resonance 
[1, 0, @, 0, H, 111- A major advantage of photoassoci- 
ation resonances is that the light intensity can be var- 
ied on short length and time scales, thus offering more 
flexible experimental control over the scattering length. 
The problem with photoassociation resonances is that 
the light induces inelastic collisions between atoms which 
lead to rapid loss of atoms. This is why photoassociation 
resonances have only rarely been used to tune the scat- 
tering length 0, Hi ■ A solution for this problem exists for 
alkali earth atoms where narrow intercombination lines 
allow for tuning of a with only moderate loss [l^, 
But this is not feasible in the large number of experi- 
ments with alkali atoms. 

In a recent experiment we explored an alternative 
scheme for controlling the scattering length with laser 
light. This scheme uses the existing coupling between an 
atom-pair state |a) ® |a) and a molecular state \g) near a 
Feshbach resonance, as illustrated in Fig. [T] By adding 
a light field that is somewhat detuned from a bound-to- 
bound transition between state \g) and an electronically 
excited molecular state |e), one can induce an ac-Stark 
shift of state \g). This results in a shift of the magnetic 
field Bres at which the Feshbach resonance occurs. If the 
magnetic field B is held close to the Feshbach resonance, 
then spatial or temporal variations of the light intensity 
affect the scattering length. This scheme also suffers from 



light-induced inelastic collisions. However, in Ref. [12 1 
we reported a two-body loss rate coefficient as small as 
K2 ^ 10^^^ cm'^/s at parameters where the real part of 
the scattering length is changed by Re(a)/abg — 1 = ±1 
with respect to its background value abg. This represents 
a reduction of the loss rate by one order of magnitude 
compared to a photoassociation resonance in ^^Rb where 
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FIG. 1: Level scheme. Pairs of atoms, each in state \a), 
are coupled to a dimer state \g) belonging to the electronic 
ground state. This coupling has a strength a and causes a 
magnetic Feshbach resonance. Laser light is applied to drive 
a bound-to-bound transition from state \g) to an electroni- 
cally excited dimer state |e) with Rabi frequency i7_R. This 
causes an ac-Stark shift (not shown) of state \g) which leads to 
a shift of the magnetic field at which the Feshbach resonance 
occurs. Typically, the same light can also drive photoassocia- 
tion from state \a) ®\a) to state |e) with coupling strength /3. 
The photoassociation is a nuisance for the scheme used here. 
The light frequency is typically somewhat detuned from both 
transitions. 



K2 ^ 10 cm^/s was reported for the same change of 
Re (a) 0,i. 

Here we show experimentally that a value as low as 
K2 ~ 10~^^ cvD? /s at Re(a)/abg —1 = 1 can be reached 
with our scheme, thus reducing the loss rate by one more 
order of magnitude compared to our previous results. 
This improvement is achieved by increasing the detun- 
ing and the intensity of the laser light, which is a com- 
monly used method for reducing incoherent rates in ex- 
periments, which rely on the ac-Stark shift. When follow- 
ing this approach, two new issues have to be addressed. 
First, the detunings are no longer small compared to the 
typical splitting between the various hyperfine and mag- 
netic substates of the excited state. Achieving a large 
detuning with respect to all excited states thus requires 
knowledge of the positions of all nearby excited states, so 
that excited-state spectroscopy has to be performed, on 
which we report in Sec. lIII Al Second, in our previous ex- 
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periment [iJl we used a model for the temporal evolution 
of the cloud size to determine Re(a) and K2- If the laser 
that drives the bound-to-bound transition is now oper- 
ated at much larger power, then it creates a noticeable 
dipole trap. Its rapid turn-on induces large-amplitude 
oscillations of the cloud size which are difficult to model. 
We therefore explore an alternative way of measuring 
Rc(a) and K2. To this end, we pin the positions of the 
atoms with a deep optical lattice which minimizes the 
effect of the additional dipole trap. In the lattice, we 
use excitation spectroscopy and a loss measurement to 
determine Re(a) and K2^ respectively. These results are 
reported in Sec. IIII Bl 

In Ref. we also studied the behavior of the system 
when the laser light is tuned close to the bound-to-boimd 
resonance. In this regime, we observed an Autler-Towncs 
doublet in the loss rate coefficient as a function of mag- 
netic field K2{B). In the same work, we presented a sys- 
tematic study of the dependence of the positions i?,.os of 
the loss resonances on laser power and detuning. In Sec. 
IlllCl of the present paper, we complement these measure- 
ments with a systematic study of the height and width 
of the resonances in this Autler-Townes doublet. These 
experimental data agree well with the theoretical model 
that we use, thus showing that the relevant physics is 
well understood. 

Before turning to the experiment, we begin in Sec. |TT] 
with a detailed discussion of the theoretical model that 
we use to describe our experiments. 



II. THEORY 

We consider a system with three internal states, as 
shown in Fig. [T] State \a) represents a single unbound 
atom. States \g) and |e) represent a dimer in the elec- 
tronic ground and excited state, respectively. The atom- 
molecule coupling between the atom-pair state \a) ® \a) 
and the dimer state \g) relevant for the Feshbach reso- 
nance is described by a coupling strength a. Laser light 
couples states \g) and |e) on a bound-to-bound transition 
with a Rabi frequency 17/}. The same laser light drives 
photoassociation from state \a) ®\a) to |e) characterized 
by a coupling strength (3. 

This model is closely related to previous studies of the 
combination of a Feshbach resonance with a photoasso- 
ciation resonance, see e.g. Refs. [H, [13, [l3| . Unhke those 
references, we are mostly interested in the Feshbach reso- 
nance and the bound-to-bound transition. The photoas- 
sociation is a nuisance in our scheme, because any useful 
change of Re(a) that it induces is inevitably accompanied 
by the loss rates that limited the photoassociation exper- 
iments in Refs. 0,i]. Luckily, we find excited states in 
our experiment for which photoassociation is negligible. 

Before turning to a quantitative model, we present a 
qualitative argument that motivates why creating a given 
change in Re (a) with our scheme causes a smaller loss 
rate than a photoassociation resonance would do. To ex- 



plain this, we must first understand the limitation of the 
photoassociation resonance. If an infinite amount of laser 
intensity were available, then the value of K2 caused by 
a photoassociation resonance could be reduced without 
modifying the change in Re(a). To this end, one would 
simply have to increase Ag along with the laser intensity, 
see Eq. ([22]) . The same improvement could be achieved if 
a photoassociation resonance with a larger transition ma- 
trix element was available because that would be equiv- 
alent to having more laser intensity. 

Our scheme relies on the fact that the typical inter- 
atomic distance is orders of magnitude larger in the 
atomic gas than within a single molecule. As a result, 
the transition matrix element is typically orders of mag- 
nitude smaller for a photoassociation resonance than for 
a bound-to-bound transition, sec Eq. Hence at a 

given laser intensity, one can detune the laser pretty far 
in our scheme and still achieve a significant light-induced 
change in Re (a), whereas a photoassociation resonance 
driven with the same laser intensity requires a smaller 
laser detuning, which results in a larger loss rate. 



A. Hamiltonian 

According to Refs. [l^, the system with the light 
off is described by the Hamiltonian 



H = 



je{a,g,e} 
1 

2 



^ / d3x*t(x)ij^.(x)vl/,(x) 



i j d3xi(i3x2^t(x^^X2)C/bg(xi2)*p(Xi,X2) 



h / d^xid^X2 



\E'J(Xi2)a(xi2)«'p(xi, X2) + H.c. 

(1) 

with the bosonic field operators ^j(x) for j e {a,g, e}, 
the single-particle Hamiltonians Hj{x.), the relative co- 
ordinate X12 = Xi — X2, the center-of-mass coordinate 
X12 = (xi -|-X2)/2, and the abbreviation >I'p(xi,X2) = 
^'a(xi)^a(x2) for the annihilation of an atom pair. Elas- 
tic two-atom collisions far away from the Feshbach reso- 
nance are described by C/bg(xi2) and the atom- molecule 
coupling relevant for the Feshbach resonance is described 
by a(xi2). 

Collisions between more than two atoms were ne- 
glected here assuming that the interatomic potentials 
are short ranged compared to the average interatomic 
distance. Furthermore, collisions involving at least one 
molecule were neglected, assuming that the density of 
molecules is low. 

a(xi2) is the Fourier transform of [l^ 



aik)^y^{g\H\k) 



(2) 



where V is the quantization volume, |k) is a two-atom 
scattering state that has an incoming plane wave with 
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wave vector k in the relative coordinate, and H is the 
Hamiltonian describing the atom-molecule coupling in 
the collision of two atoms. Its position representation 
has the asymptotic form 



ikxi 



(Xl2|k) 



/(k)e^'=--/xi2 



Xl2 



oo, 



(3) 



where /(k) is the scattering amplitude. Combination of 
the last two equations shows that a(k) is independent of 
V. 

If all incoming particles in the gas are slow enough, 
then we can replace Q!(k) by a constant a = linik-^o a(k). 
In the position representation, this corresponds to a con- 
tact potential in the relative coordinate. We use an anal- 
ogous approximation for Uhg and obtain H = J d^x7Y(x) 
with the Hamiltonian density [3| 



H(x) 



E 

je{a.g,e} 



^ *](x)iJ,(x)#,(x) 



h [*1,(x)a#p(x) +H.C. 



(4) 



Note that Uhg — A-nh'^ahg/'m is related to the background 
scattering length abg and the atomic mass m. 

We now extend this model to include effects caused 
by the light, which has an electric field of the form 
E = — Ef) cos{uj i,t) with amplitude Eq and angular fre- 
quency ojl- In analogy to Eq. ^ we obtain the following 
additional terms for the Hamiltonian density 



Hl(x) = hcos{ujLt) *l(x)r!fi*g(x)+H 



+ 2ncos{uLt) *1:(x)/3#p(x) -hH.c. . (5) 

The Rabi frequency and the coupling strength (3 de- 
scribe the bound-to-bound transition and the photoasso- 
ciation, respectively, is related to a matrix element 
of the matter-light interaction term in the Hamiltonian, 
which in the electric dipole approximation reads Hit) = 
-dE{t). This yields hnRCOs{ujLt) = -{e\d\g)E{t). We 
abbreviate c?eg = and obtain the well-known rela- 

tion 



End 



0"eg 



For P we obtain in analogy to Eq. ^ 



(6) 



(7) 



Finally, we specify the single-particle Hamiltonians 
which consist of a kinetic part and an internal part E™^ 



i7,(x) 



2m o 



E" 



(8) 



for j € {a,g,e}. Here rrij = m(2 — 6aj) is the mass of a 
particle in state j and Sij is the Kronecker symbol. The 



internal energy of each state depends nonlinearly on the 
magnetic field B. Near the pole of the unshifted Feshbach 
resonance Bpoio we approximate this dependence as lin- 
ear and obtain 



E: 



int 



~^iJ{B-B^ 



(9) 



where /Xj is a magnetic dipole moment. At B = Bpoio, 
the internal states |a) and \g) are degenerate whereas the 
internal state |e) has an energy offset fiujeg- 



B. Mean-Field Model 

We assume that the population in each state is Bose 
condensed and can be described by a mean field ijjj (x) = 
(\E'j(x)). In addition, we assume that the system is homo- 
geneous. We take the expectation value of the Heisenberg 
equation of motion ifid^ j / dt = j, H] and approximate 
the field operators as uncorrelated, so that the expecta- 
tion values factorize [Tsj . Physics beyond this approxi- 
mation is discussed in Refs. [itI. [Toj. We obtain 



/ 



pint jj 

+2a*V:V's + 4/3* cos{ujLt) (10a) 
i-^-ipg = ai;l + ^i:g + n*j^ipeCOs{uJLt) (10b) 



(2/?V2 ^ ^^^^^ coa{LOLt) 



E. 



int 



-V'e.(lOc) 



We move to an interaction picture by replacing ip,, — > 
^^^itAJLt g^j^jj perform a rotating- wave approximation by 
neglecting coefficients rotating as g^^*'^^*. We then 
move to another interaction picture by replacing ijjj — s- 
4>j exp ((2 — 6ni)iE™*^t/h) for j G {a, g, e}. Hence (see also 
Ref. [12) 



dt'^'' h 



IVal + 2a>:Vg + 2/3», (11a) 



A. 



where we abbreviated 
1 



Ai 



f^ag{B — -Bpolc) 



-Al + —Hae{B — Spolc) 

n 



eg 



(lib) 

7e ) Ve, (lie) 



(12a) 

(12b) 
(12c) 



with ^ag = 2/ia — ^'^d ^ae = ^^a— IJ-e- In Eq. (jllcP WC 

included an ad hoc decay rate 7e that represents spon- 
taneous radiative decay from state |e) into states that 
are not included in the model, similar to R ef. l20l| . A 
model similar to Eqs. pT|) was used in Ref. [Tal to ex- 
plain enhanced photoassociation loss rates near a mag- 
netic Feshbach resonance 
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Note that the typical interatomic distance is orders of 
magnitude larger in the atomic gas than within a single 
molecule. For a typical excited state, this results in 



(13) 



C. Adiabatic Elimination and Scattering Length 

We assume that all the population is initially prepared 
in state \a) and that the populations in states \g) and |e) 
will remain small at all times so that theycan be elim- 
inated adiabatically, similar to Refs. [IB, H^]. This is a 
good approximation, e.g., if the angular frequencies a^a 
and (3^a are both small compared to fi^j and 7e or com- 
pared to Ae and Ag. This condition is always satisfied 
in the low-density limit, but for a very broad Feshbach 
resonance it might be difhcult to reach this regime ex- 
perimentally. 

The adiabatic elimination is achieved by formally set- 
ting (d/dt)ipg = (d/dt)ipe = 0. This is used to eliminate 
ipg and ipe from the equations. We obtain 



A-nTia 



m 



(14) 



with the complex-valued scattering length 
m 

|a|2(Ae-^7e/2)-Re(a*^!^J/3)^ 



a 



X 



^A„ 



(Ae - lJe/2)Ag - |17fl/2|^ 



(15) 



The term Re(a*ri^/3) represents interference between the 
two possible ways to go from state \a) to state |e), either 
directly or indirectly through state \g). 

The real part of the scattering length is responsible for 
the mean- field energy [2lj. We assume that Obg is real 
and obtain 



Re(a 



2TTh (AgAe - |r!fl/2|2 



(As7e/2) 



A„ A,A„ 



a 



'A,(7e/2)^ 



(16) 



The imaginary part of the scatterin g le ngth gives rise 
to two-body loss with a rate equation 



dn 
'dt 



2„(2) 



im(a), 

m 



(17a) 
(17b) 



where n = j-^ap is the atomic density, K2 is the two-body 
loss coefficient for a Bose-Einstein condensate (BEC), 
and 5*-^' is the pair correlation function at zero relative 
distance. For a BEC with N atoms g'^^ = 1 - 1/iV. 
Insertion of Eq. (fT5|) yields 



270 



(AgAe - |f^fl/2|2)2 + iAgjj2y 



(18) 
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FIG. 2: (Color online) Predictions for Re(a)/a6g and A'2 as 
a function of B from Eqs. (|16|l and (|18p . For all curves, we 
choose h-ye/fiae = 2 G, /3 = 0, and AB = 0.2 G; with AB 
defined in Eq. (|20|l . The dotted lines (red) are for resonant 
laser light Al = and TilOnf /jelJ-ag ~ 100 G. They show 
two resonances that are symmetrically split around Bpoic. 
These resonances represent an Antler- Townes doublet. The 
solid lines (blue) are for large laser detuning AL/7e = 5 and 
?i|r2_Rp/7e^a9 = 100 G. They also show two resonances, but 
their heights, widths, and distances from -Bpoic are quite dif- 
ferent. The dashed line (black) is a reference without any 
light Qr = 0. 



All terms in the numerator are oc Eq, so that the relative 
importance of the terms is independent of laser inten- 
sity. If B is held near the unshifted Feshbach resonance, 
then Ag is small and all terms containing f3 in Eq. 
become negligible. Hence if one considers measurements 
of K2 performed fairly close to the Feshbach resonance, 
then photoassociation is negligible and it is impossible 
to extract the value of f3 only from such measurements 
(unless Qn vanishes). 

The interference term in Eq. (jlSp leads to a correspond- 
ing interference term in Eq. (|18p. Note that the minimum 
of K2 observed in Ref. [ij] is a result of destructive in- 
terference due to this term [Tst . 

Figure [5] shows predictions for Re(a)/abg and K2 as a 
function of B. For large \flR\, one can clearly see two res- 
onances in K2 each of which is approximately Lorentzian. 
Each of these resonances is accompanied by a dispersive 
feature in Re(a)/af,g. 

We show now that our model reproduces known results 
from the literature in the special cases of a pure Feshbach 
resonance or a pure photoassociation resonance. For a 
pure Feshbach resonance {(3 = ^In = 0) Eq. yields 
the familiar result [l| 



1 



AB 



B-B, 



pole 



(19) 
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with the width of the Feshbach resonance 



AB 



(20) 



For a pure photoassociation resonance [a = CIr = 0) Eq. 
(Uni) yields 



a — Qhg 



2^?% Ae - i7e/2' 



(21) 



which is a Breit-Wigner form [22] as a function of or 
B. The real and imaginary parts are 



Rc(a) = flbg - T^A, 

l/^|2 

K2 = 27, 



27r?i"^A2 + (7,/2)2 ^22^) 



^A2 + (7,/2)2' 



(22b) 



which is identical to Eq. (10) in Ref. @ in the limit 

Estim Egpon- 

A quantitative comparison of E g. (1181) with our ex- 
perimental data (here and in Ref. [l3]) shows that /? is 
negligible for most of the excited states |e) that we use. 
We therefore set /3 = for the rest of the theory section. 



D. Large Detuning 

A good part of our experiments is performed in the 
limit of large laser detuning where \Al\ ^ |/iae(-B — 
■Bpoio)/?i| and |Ai| > and with /3 = 0. In Eq. pS]) 
we can approximate l/(Ae — i7e/2) ~ — l/A^ — i7e/2A| 
and obtain a Breit-Wigner form 



Obg 1 



AB 



B — Bjc 



iW/2 



(23) 



with AB from Eq. (j20p . The real and imaginary parts 
are a dispersive line shape and a Lorentzian, respectively, 



Re(a) = atg ( 1 - ) (24a) 



K2 



{B - B,,,)^ + wyA, 

^max 

^2 (24b) 



1 + AiB - B,es)yW^ 



The resonance position B^-^s, the maximum loss rate co- 
efficient K2^^^, and the full width at half maximum W 
of the Lorentzian are given by 



r^max 
J<2 



w 

-^rcs ^ -^polc 



_h_ 8|ap 



Ma, 4Ai^^ 
Mag 4Ai 



(25a) 
(25b) 
(25c) 



The far-detuned bound-to-bound coupling yields the 
well-known ac-Stark shift of state \g) and this shifts B^^s- 
As in most applications of ac-Stark shifts, we wish to 
achieve a certain value of jfiijp/Ai and at the same time 
keep the rates for incoherent processes as low as possible. 
Hence, it is advantageous to increase the detuning and 
power of the laser in a way that keeps \flji\'^/AL constant. 
This yields ^ and ^2(5) ^ {ATTh\a\^ / ^iag)SiB - 
Bics), where d denotes the Dirac delta function. For any 
given value of i? ^ i?ios one can thus decrease K2{B) by 
increasing the detuning and the laser power sufRciently 
far. 

In general, it is possible that several excited states con- 
tribute noticeably to a. Our model is easily adapted to 
this situation by introducing a separate version of Eq. 
(|llcp for each excited state and by including sums over 
the excited states in Eqs. (jllap and (jllbp . In the limit 
of large laser detuning and with /3 = for each excited 
state, Eqs. (|23p - (|25a[) remain unchanged and a sum over 
the excited states appears on the right hand side of Eqs. 
(I25bl) and (1^ . 



E. Autler-Townes Model for Weak Damping 

More insight into the physics of the problem can be 
gained from an Autler-Townes model [2^, [IJ]. In ad- 
dition, analytic expressions for the position, height and 
width of the resonances in K2{B) can be derived. 

This approach is based on the assumption that the 
dominant frequencies in the problem are Qr and/or 
(Ag — Ag). In this case, one can first diagonalize the 
driven two-level system spanned by \g) and |e) and sub- 
sequently treat the coupling to state \a) as a weak probe. 

For the first step, we diagonalize the two-level system 
spanned by \g) and |e), setting a = /3 = 7e = 0. We 
assume without loss of generality that the relative phase 
between states \g) and |e) is chosen such that flu is real. 
This yields energy eigenvalues and eigenvectors 



E± = -(A, + Ag±r!eff) (26a) 

1+) = cos^|e)+sin^|5) (26b) 

h) = -sin-|e) -f-cos-lg), (26c) 

where the effective Rabi angular frequency Qcff and the 
mixing angle ■(? are real-valued and must satisfy the im- 
plicit equations 



fJoffCOSl? : 

f2cff sin 'd = Q 



(27a) 
(27b) 



This determines a unique value of i? modulo 27r and it 
yields 



a 



off 



(Ae 



Ag)2. 



(28) 
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For the second step, we rewrite the mean-field model 
(fTTj) in the new basis and obtain with (3 = 



/ 

dt 


h ' 




fC^T/.-) (29a) 


d 


= c+4i^ 




- ^7mixV'- (29b) 








- ^7_^ (29c) 


with C-f 


= a sm -J , 


C_ = a COS ^, 7+ 


= 7eC0S^ |, 7_ = 


7eSin^ |. 


, and 7niix 


= — 7eSm^cos^ 


. So far wc only 



rotated the basis and the model is still exact. We now 
approximate the loss as being diagonal in the states |+) 
and |— ) by setting 7mix = 0. Adiabatic elimination of 
the populations in states |+) and |— ) then yields 



abg + a+ + a- 



(30) 



where the states |+) and |— ) each contribute a Breit- 
Wigner form as a function of E± 



a± 



\C±\ 



2Trh E±/h-ij±/2' 



(31) 



Thus, the states |+) and |— ) each cause a single res- 
onance and their contributions to a are simply added. 
This corresponds to the intuitive understanding of an 
Autlcr-Townes doublet [24j . 

The above approximation 7niix = is self-consistent if 
the system is close to one resonance and the resonances 
are well separated (7e ^ ^^off), because in this case the 
states |±) have very different populations, so that a pos- 
sible coherence between these populations has little effect 
and 7,nix is negligible. 



F. Properties of the Autler-Townes Resonances 



We note that E±, 7±, and C± in Eq. ([3T|) are gener- 
ally nonlinear functions of B. But if the resonances are 
narrow and well separated, then will vary only little 
within the width of a resonance and we can approximate 
7± and C± as constant across a resonance. According 
to Eqs. (|17b[) and (|3ip . maxima of K2{B) will then oc- 
cur at E± = 0. Combination with Eqs. (|12ap . (|12b|) . 
P6ap . and (|28p yields the magnetic fields i?res at which 
the resonances occur 



Bics -^polc ^ 



The condition E± = combined with Eqs. (jlTbp and 
3T|) yields the maximum of the loss rate coefficient 



A— =^rtan^ 
7e V 2 



±2 



(33) 



Insertion of cos^ | = 5(1 + cost?), sin^ f ^ 5(1 ^ cosi?), 
Eqs. (P5a|) . ((77a|) . and E± = yields 

lap Ae 



^max 



(34) 



7e Ag 

Insertion of Eqs. (fna|) . (fT2bl) . (15^ . and B = B,.cs yields 



^ 7e A'ag 



2A, 



(35) 



The width of the resonance can also be calculated. For 
this, we recall that the resonance is narrow and approx- 
imate E± as linear in B around i?rcs- As a result, Eq. 
([3T|) becomes a Breit-Wigner form as a function of B. 
Correspondingly K2{B) becomes a Lorentzian as in Eq. 
(I24bl) . 

In this linear approximation E± ~ e{B — -Brcs) with 
e = dE±/dB\B,,^. Eqs. p^ . p2bl) . (pSa|l and dMl) yield 



e = 2 (Mae + Mag ± (/^ae " A'ag) 



(36) 



According to Eqs. (jl7bp and the fuU width at half 
maximum of the Lorentzian K2{B) isW = ?i7±/e. Using 
Eq. (|27ap we obtain 



n-ie 2e f § 

w^ = T±^='^'^^+^- r"2 



±2 



(37) 



The term (tan |) ' can be evaluated as in Eq. (|33p . 
yielding 



A, 



A 



2 ^ iia^f^2 



Mag 



(38) 



We will compare these results with experimental data in 

Sec. Unci 



III. EXPERIMENT 

Our experiments are performed with ®^Rb using a 
Feshbach resonance that is characterized by the param- 
eters m, m m m ^poio = 1007.4 c, as = 0.21 c, 

^lal^-Kh. = 1.02 MHz/G, fXag/^TTh = 3.8 MHz/G, and 
flbg = 100. 5ao, where ag is the Bohr radius. With these 
parameters Eq. ^ yields |Q|/27r = 1.8 mHz cm^/^. The 
peak density of the BEC is typically n = \4'a\'^ ~ 2 x 10^'' 
cm"3. The quantity |%/8aVa| = 2tt x 0.07 MHz can be 
regarded as a Rabi angular frequency. For typical param- 
eters of our experiment, this value is small compared to 
7e and \^r\-, so that the adiabatic elimination performed 
in Sec. Ill Cl is justified. To set the scale for K2 in Fig. [51 
we note that Uhg/h = 4.9 x 10~^^ cm^/s. 




2 3 4 5 6 7 8 

laser frequency (GHz) - 382040 

FIG. 3: (Color online) Excited-state spectroscopy for *'Rb. (a) Photoassociation spectrum taken at B = 1000.0 G. (b) 
Bound-to-bound spectrum taken at B = 1000.0 G. Arrows indicate the resonances characterized in Tab. IT] 



A. Excited- State Spectroscopy 

The starting point for all experiments in this paper is 
an essentially pure BEC of ^^Rb atoms in the hyperfine 
state \F,mp) = These atoms are held in a crossed- 

beam optical dipolc trap with both beams operated at 
1064 nm and with trap frequencies of (ujx,iOy,iOz)l'2.i: ~ 
(74, 33, 33) Hz. Gravity acts along the x axis. The mag- 
netic field B points along the z axis and is held several 
gauss away from the Feshbach resonance. 

For the excited-state spectroscopy, we use two tech- 
niques which complement each other. The first technique 
is ordinary photoassociation spectroscopy. For this, we 
simply illuminate the BEC with photoassociation light. 
The light intensity is slowly increased within 80 ms to 
a final power of ^ 10 mW and held there for 100 ms. 
The slow increase is needed to avoid large-amplitude os- 
cillations of the cloud shape. Next, the photoassociation 
light, B, and the dipole trap are switched off simultane- 
ously. Finally, the remaining atom number is extracted 
from a time-of-flight image. 

The photoassociation light is implemented as a 
traveling-wave laser beam with a waist (1/e^ radius of 
intensity) oi w ~ 0.17 mm and a wavelength of ^ 784.7 
nm. In order to address as many excited states as possible 
in the spectroscopy measurements, we let the photoas- 
sociation beam propagate along the x axis and choose 
a specific linear polarization which corresponds to 1/3 
of the intensity in each of the polarizations tt, (t+, and 
(7^. As the population is initially in the atomic state |a), 
the loss signal is typically dominated by photoassociation 
processes so that the technique is particularly sensitive 
to excited states that have a large value of the photoas- 
sociation coupling strength /3. 

Figure a) shows a photoassociation spectrum sX B ~ 
1000.0 G. We extended this photoassociation scan down 
to 382,037.3 GHz but did not find any further loss reso- 



nances. We determined the corresponding zero-field fre- 
quencies by performing a similar measurement at _B ~ 0. 
We took data between 382,034 GHz and 382,051 GHz 
and found photoassociation loss resonances in the range 
between 382,041.8 GHz and 382,044.5 GHz. A compari- 
son with photoassociation data from the Heinzen group 
shows that the excited states involved are the hyperfine 
and magnetic substates of the vibrational state v = 120 
in the attractive Ig potential that is adiabatically con- 



nected to the P- 



3/2 



Si/a threshold [2S 



Since our technique for shifting a Feshbach resonance 
with laser light relies on a bound-to-bound transition, not 
on photoassociation, we developed a second spectroscopy 
technique that is particularly sensitive to excited states 
with a large value of VLr- The basic idea is to first use 
the Feshbach resonance to associate molecules into state 
\g) and then illuminate them with light that resonantly 
drives bound-to-bound transitions. We call this the 
bound-to-bound light and employ the same laser beam 
previously used for the photoassociation spectroscopy. 

In order to avoid loss of particles due to inelastic colli- 
sions between molecules, the atoms must be loaded into 
a deep optical lattice before associating the molecules 
[S^l- The lattice has a light wavelength of 830.440 nm 



TABLE I: Parameters of four selected bound-to-bound reso- 
nances, e is the elementary charge, ao the Bohr radius. /3 is 
negligible for all these resonances. 



Polarization 


LJeg/27r 




7e/27r 


^ae/27rft 




(MHz) 




(MHz) 


(MHz/G) 


cr~ 


382,045,759.4(3) 


0.24(5) 


4.4(5) 


2.2(1) 


a+ 


382,045,818.2(3) 


0.29(5) 


4.7(5) 


1.7(1) 


TT 


382,046,942.8(3) 


0.28(5) 


4.7(5) 


2.6(1) 


TT 


382,047,581.8(3) 


0.18(5) 


5.3(5) 


2.7(1) 
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and a depth of Vq ~ 20Er, where Er is the atomic recoil 
energy. As in Ref. [3l| we prepare an atomic Mott insu- 
lator, which contains exactly two atoms at each lattice 
site in the central region of the lattice. This core is sur- 
rounded by a shell of sites that contain exactly one atom 
each. After loading the lattice, the laser power of one of 
the dipole trapping beams is ramped to zero, as in Ref. 

Next, we ramp the magnetic field slowly downward 
across the Fcshbach resonance so that molecules are as- 
sociated [SLii^] in state \g) at sites that contain exactly 
two atoms. Sites containing three or more atoms might 
exist due to imperfect state preparation. Such sites are 
emptied by inelastic collisions as soon as molecules begin 
to form. Sites containing one atom are unaffected by the 
ramp. 

Next, the bound-to-bound light is turned on for 0.2 ms 
at a power of ~ 0.1 fiW. This light has the same linear 
polarization as for the photoassociation spectroscopy. If 
a molecule in state \g) is excited on a bound-to-bound 
transition, then it is likely to undergo spontaneous ra- 
diative decay into a different internal state. After turn- 
ing off the bound-to-bound light, the magnetic field is 
ramped back across the Feshbach resonance to dissociate 
the molecules that remained in state \g). Subsequently, 
the optical dipole trap at 1064 nm is turned back on, 
the lattice depth is slowly ramped to zero, the cloud is 
released, B is switched off, and the atom number is deter- 
mined from a time-of-flight image. Molecules that were 
excited by the laser are not dissociated and thus not de- 
tected. 

A bound-to-bound spectrum measured at B = 1000.0 
G is shown in Fig. [3jb). Comparison with the photoas- 
sociation spectrum in part (a) shows that many of the 
excited states are visible with both techniques. But iden- 
tifying promising candidates with large is not easily 
possible from part (a). The light frequency calibration 
for both spectra has a precision of ~ 30 MHz and can 
fluctuate within a single scan. This causes deviations 
in the resonance positions between the photoassociation 
spectrum and the bound-to-bound spectrum. In addi- 
tion, the states |a) and \g) are degenerate at B = 1007.4 
G. At B = 1000.0 G, the internal energies of a molecule 
in state \g) and a pair of atoms in state \a) differ by 
2TTh X 20 MHz [131. In Fig. [3] this yields a 20-MHz shift 
of all bound-to-bound resonances with respect to the pho- 
toassociation resonances. 

Very different values of the light intensity and the illu- 
mination time were used when recording the two spectra. 
This tremendous difference in the sensitivity of the two 
methods is a result of Eq. . 

In Ref. [l3l we developed a method to determine all 
the parameters of a bound-to-bound resonance. We now 
apply this method to four reasonably strong bound-to- 
bound resonances which are fairly close to the high- 
frequency end of the spectrum in Fig. [31 Results are listed 
in Tab. HI For these measurements, the laser producing 
the bound-to-bound light was beat-locked to a frequency 



comb, resulting in a much better precision of the fre- 
quency calibration. The polarization of each resonance 
was determined from a series of measurements in which 
the bound-to-bound light had only one of the polariza- 
tions TT, CT"*", or a~ . The latter two polarizations were 
implemented with the bound-to-bound beam propagat- 
ing along the z axis. Each resonance in Tab. U responded 
to only one of these polarizations. 

Our choice of the light wavelength for the experiments 
in the following Sec. IIII Bl is based on the spectra ob- 
tained here. In order to achieve a large detuning from 
all excited states, the light must be detuned either to the 
left or to the right of the complete spectrum in Fig. [31 
The three outermost resonances at the low-frequency end 
of the spectrum in Fig. [31 have fairly strong photoassoci- 
ation loss features, but show hardly any bound-to-bound 
features, which is unfortunate. The high-frequency end 
of the spectrum looks more promising. Wc therefore per- 
form all the following experiments blue detuned from the 
high-frequency end of the spectrum in Fig. [31 According 
to Tab. m the two strongest bound-to-bound resonances 
near this end of the spectrum both respond to tt polarized 
light. Hence, we choose n polarization for the bound-to- 
bound light in all the following experiments. 

We note as a side remark, that our model, unlike Ref. 
p^ . assumes the existence of a direct bound-to-bound 
coupling term O^. In Ref. [isj an indirect bound-to- 
bound coupling is constructed by invoking virtual tran- 
sitions into the continuum of excited atom-pair states 
above threshold. This implies that the bound-to-bound 
coupling should be proportional to the photoassociation 
coupling P [l^ . Our spectroscopy data in Fig. [31 do not 



support this prediction of Ref. [15[. There are strong 
photoassociation resonance that have hardly any bound- 
to-bound coupling. This indicates that the indirect cou- 
pling discussed in Ref. [ll] is neghgible [ssj . 



B. Shifting the Feshbach Resonance with Light 

Wc now use the spectroscopic information gathered 
above to shift the Feshbach resonance with far-detuned 
light. As discussed in the introduction of this paper, 
we minimize the effect of the dipole trap created by the 
bound-to-bound light by working in a deep optical lat- 
tice. 

In order to measure Re(a), we first load the atoms into 
the lattice as described in Sec. IIII Al We then use exci- 
tation spectroscopy [s^ in the lattice, i.e., we modulate 
the power of one rctro-reflcctcd lattice beam sinusoidally 
as a function of time around an average lattice depth of 
Vo ^ IbE,.. The modulation amplitude is ~ iE,.. The 
modulation lasts for 10 or 20 ms. During the modulation, 
the atoms are illuminated with the bound-to-bound light 
and B is held at a specific value close to the Feshbach 
resonance. The bound-to-bound light is on for a long 
enough time that sites containing two or more atoms are 
essentially emptied by light-induced inelastic collisions. 
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The signal in the excitation spectrum that is sensitive 
to the modulation frequency thus stems from sites that 
were initially populated by one atom. For certain modu- 
lation frequencies, tunneling of an atom between two such 
sites is resonantly enhanced. This leads to a frequency- 
dependent loss of atoms and of atomic phase coherence. 
At the end of the modulation, wc switch the bound-to- 
bound light off and simultaneously jump B back to a 
value several gauss away from the Feshbach resonance. 
Next, the dipole trap at 1064 nm is turned back on and 
the lattice depth is slowly reduced to Vq ^ 6Er, where 
the gas is superfluid, thus restoring phase coherence be- 
tween neighboring lattice sites. Finally, the dipole trap, 
B, and the lattice are simultaneously switched off. The 
time-of-flight image shows satellite peaks due to the re- 
stored phase coherence. 

The visibility [IHl of the satellite peaks displays a min- 
imum at a modulation frequency where tunneling pro- 
cesses between two initially singly occupied sites are res- 
onant. This minimum is located at a frequency / = 
Ke{U)/2T:h with the on-site interaction matrix-element 
U ~ g J (i'^a;|u'(x)|^, where g = Airh^a/m and w is a 
tight-binding Wannier function. The measurement of / 
thus yields Re(a). 

A sequence of such measurements for various values 
of B yields Fig. [UJa). For parameters where Re(a) is 
reduced drastically, the system becomes superfluid and 
the peak in the excitation spectrum is smeared out so 
much that its center cannot be determined any more. 
Hence, this method is not applicable in this regime. For 
comparison, the figure also shows Re(a) measured with 
the same method, but in the absence of bound-to-boimd 
light. Clearly, the position B^cs of the Feshbach reso- 
nance is shifted by ~ —0.35 G due to the presence of the 
light. 

We now turn to the question, how large a loss-rate co- 
efficient K2 is associated with this shift. In order to de- 
termine K2, we load the atoms into the lattice and asso- 
ciate molecules as described in Sec. IIII Al Right after as- 
sociation, we dissociate the molecules. This association- 
dissociation sequence serves the purpose of emptying all 
sites that contain three or more atoms, which will be- 
come important below. The lattice depth is Vq ~ 20Er 
so that tunneling is negligible. We then switch on the 
bound-to-bound light and simultaneously jump B to a 
value close to the Feshbach resonance. These conditions 
are maintained for a variable hold time. During this hold 
time, the bound-to-bound light causes rapid loss of atoms 
in doubly occupied sites. Next, B is switched to a value 
several gauss away from the Feshbach resonance and the 
bound-to-bound light is switched off. The dipole trap 
light at 1064 nm is turned back on and the lattice depth 
is slowly lowered to zero. Finally, the cloud is released 
and B is switched off. The remaining number of atoms 
is extracted from a time-of-flight image. 

The two-body loss during the hold time is described by 
the master equation of Ref. [1^. We consider a lattice 
site initially occupied by exactly two atoms and we ne- 
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FIG. 4: (Color online) Shifting the Feshbach resonance with 
laser light, (a) Elastic and (b) inelastic two-body scattering 
properties are shown as a function of magnetic field B. Ex- 
perimental data in the presence (•) and absence (o) of the 
light are compared. The light power is 11.2 mW and the fre- 
quency is L0L/2n = 382, 048, 158 MHz, which is 576 MHz blue 
detuned from the nearest bound-to-bound transition. The 
solid lines in (a) and (b) show fits of Eqs. (|19|l and (|24b|) . re- 
spectively, to the data. The Feshbach resonance is shifted by 
~ -0.35 G. At B = 1006.91 G, we measure Re(a) /obg - 1 ~ 1 
and 7^2 ~ 1 X 10"^^ cm^/s. 



gleet tunneling between sites. The master equation then 
yields the density matrix p = p\2) (2| -|-(1— p)|0) (0|, where 
p = exp(— Ft) is the probability that a decay at this site 
occurred and \n) denotes a Fock state with n atoms. The 
parameter F is given by [H, [sj 



r = K2 d^x\wix)\^ 



(39) 



The decay of the total atom number N in the experiment 
is obtained by taking the sum over a large number of 
isolated lattice sites, yielding 



iv(i) = 7Vi + iV2cxp(-ri), 



(40) 



where A'^i and N2 arc the initial atom numbers on singly 
and doubly occupied sites, respectively. It is crucial that 
there are no sites with three or more atoms, because they 
would give rise to an additional term that would decay 
more rapidly, thus making it more difficult to extract F 
from the measured N{t). 

Wc measured N{t) at a fixed hold time t = 2.1 ms for 
various values of i? and used Eqs. ([55]) and (I^O]) to extract 
K2{B). Resuhs are shown in Fig.g^b). At B = 1006.91 
G, we measure Re(a)/abg — 1^1 and K2 ^ ^ ^ 10^^^ 
cm'^/s, which is one order of magnitude lower than our 
previously published result of Ref. [l3| and two orders of 
magnitude lower than the corresponding result reported 
for photoassociation resonances in ^^Rb 0, [l] ■ 
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We use the results of Sec. IIIDI to calculate the theoret- 
ical expectations from the sum of the two tt resonances 
in Tab. [H We thus expect Bpes — Bpoic ~ —0.24 G and 
K2 = 3x 10"i3 cm^/s at \B - B^csl = AS which in the 
model corresponds to Re(a)/abg — 1 = ±1- Both exper- 
imentally observed values are somewhat larger than the 
expectation. This might be due to contributions from 
other bound-to-bound and photoassociation resonances 
that we did not include in this estimate. 

We tried to reduce K2 even further by setting the 
power of the bound-to-bound light to 66 mW and its 
frequency to 382,050,911 MHz which corresponds to a 
detuning of Al/2tt = 3.33 GHz from the nearest bound- 
to-bound resonance. The expected and observed shifts 
were —0.35 G and ~ —0.65 G, respectively. But here 
we observed K2 3 x 10"^^ cm^/s at \B - B^^s] = AS, 
which is much worse than the expectation K2 ~ Ix 10~^^ 
cm'^/s. As the detuning is much larger than in Fig. HI 
other bound-to-bound and photoassociation resonances 
contribute even more strongly to the signal, which might 
explain the increased deviation between the observed val- 
ues and the estimates based on the two tt resonances of 
Tab.U 



C. Autler-Townes Doublet 

Finally, we compare the theoretical results for the 
width and height of the Autler-Townes resonances in 
K2{B) from Scc. lIIFI with experimental results. We mea- 
sured Autler-Townes doublets in K2{B), as shown in Fig. 
^h). We fit a Lorentzian (j24b[) to each of the two peaks 
in the experimental data. The best-fit values for K^^'^ 
and W are shown in Fig. [5] The experimental procedure 
and parameters are identical to Fig. 4 of Ref. [l^; see 
this reference for details. 

The experiment is performed in a regime where {^IrI 3> 
7e so that the Autler-Townes model of Sec. IIIFI is ex- 
pected to be a good approximation. The dotted lines 
show the corresponding predictions ([55]) and ([55]) . They 
agree well with the experimental data. 

For comparison, we numerically determined the peaks 
in K2{B) from the full model with the parameters 
of Tab. |T1 The corresponding maximum values K^^^ are 
shown as solid lines in Fig. [^a). As the full model 
does not predict Lorentzian lines, a direct comparison 
with the width W is not straightforward. We decide to 
use the second derivative of K2{B) at the maximum for 
a comparison. The solid lines in Fig. EJb) therefore show 
the values of 

-(-^^)""^ 

at the peaks calculated from the full model ([18]). If the 
peaks in the model were Lorentzian, this would yield the 
width W. The solid lines also agree well with the exper- 
imental data and with the dotted lines. 
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FIG. 5: (Color online) Systematic study of the loss reso- 
nances. K2{B) was measured for certain values of the laser 
frequency at a fixed laser power of 0.47 mW. (a) The max- 
imum is'™'"^ and (b) the width W were determined from a 
fit to Eq. ()24b|) . The experimental data for the resonances 
that occur at the lower (o) and higher (•) value of B both 
agree well with the predictions of the full model Eq. psp 
(solid lines) which is well approximated by Eqs. (|35|) and (|38|) 
(dotted lines). 

IV. CONCLUSION 

To summarize, we improved our recently developed 
scheme for shifting a magnetic Feshbach resonance with 
laser light by exploring the regime of even larger detuning 
and laser power. We demonstrated that the light-induced 
loss rate can be reduced by one order of magnitude com- 
pared to our pervious work [l^ . The measurements re- 
quired excited-state spectroscopy and an optical lattice. 
We also presented a detailed discussion of a model that 
describes our experimental data. 
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